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ABSTRACT 

Aims. Today, large ground-based instruments, like VISIR on the VLT, providing diffraction-limited (~0.3 arcsec) images 
in the mid-infrared where strong PAH features appear enable us to see the flaring structure of the disks around Herbig 
Ae stars. Although great progress has been made in modelling the disk with radiative transfer models able to reproduce 
the spectral energy distribution (SED) of Herbig Ae stars, the constraints brought by images have not been yet fully 
exploited. Here, we are interested in checking if these new observational imaging constraints can be accounted for by 
predictions based on existing models of passive centrally irradiated hydrostatic disks made to fit the SEDs of the Herbig 
Ae stars. 

Methods. The images taken by VISIR in the 8.6 and 11.3 /im aromatic features reveal a large flaring disk around 
HD97048 inclined to the line of sight. In order to analyse the spatial distribution of these data, we use a disk model 
which includes the most up to date understanding of disk structure and physics around Herbig Ae stars with grains in 
thermal equilibrium in addition to transiently-heated PAHs. 

Results. We compare the observed spatial distribution of the PAH emission feature and the adjacent continuum emission 
with predictions based on existing full disk models. Both SED and spatial distribution are in very good agreement with 
the model predictions for common disk parameters. 

Conclusions. We take the general agreement between observations and predictions as a strong support for the physical 
pictures underlying our flared disk model. 
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1. Introduction 

Herbig Ae/Be stars are thought to be the intermedi- 
ate mass (~ 2-10 M w ) counterparts of T Tauri stars 
(jWaters fc Waelkensl [1998). They have a strong infrared 
(IR) excess coming from warm circumstellar dust. Similar 
to T Tauri stars, this dust is believed, due to the interpre- 
tation of spectral energy distri butions (SEDs), to reside in 
a flared/flat circumstellar disk (jMeeus et al.ll2001h . 

PAHs (Polycyclic Aromatic Hydrocarbons) undergo 
transient heating : they do not reach thermal equilibrium 
with the radiation field but absorb individual photons, ex- 
periencing a rapid increase in temperature, slowly cool, thus 
re-radiating the absorbed energy in the infrared (IR) . This 
radiation allows to see much further in the disk since these 
grains can reach high temperatures far away from the star. 
As a result, PAHs emission can be use d to probe the exter - 
nal region of disks around HAe stars (jLagage et al.l 12006). 
Another point of interest is that, if present, PAHs can con- 
stitute an important source of opacity and are likely to play 
a key role in the thermal budget and chemistry of the gas, 
as they do in the interstellar medium. The thermal coupling 
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between PAHs and gas via the photoelectric effect will, for 
example, determine the gas temperature in the upper layer 
of the disk where the gas and dust temperatures are not well 
coupled (Jonkheid et al. 2004; Kamp & Dullemond 2004). 
Furthermore, PAHs are also good tracers of the presence of 
very small particles in the surface layers of disks and their 
emission can tell if and where very small particles survive 
settling and coagulation processes that cause the majority 
of the original grain population to grow to larger sizes in 
the same objects. 

We have started a program of imaging nearby HAe stars 
with VISIR, the new mid-infrared instrument attached to 
the third 8-m unit of ESO's Very Large Telescopes lo- 
cated at Cerro Paranal, Chile. One of the first targets was 
HD 97048 which has a strong IR excess, characteristic of a 
flared disk, and strong PAH features. 

In a previous work, we were able for the first time to 
constrain the geome t ry of the disk, t hanks to PAH emission 
(jLagage et al.l 120061 ) . lLagage et al.l (2006) found there is a 
large flaring disk around HD97048 extending at least up 
to 370 AU, vertically optically thick at the observed wave- 
length (8.6 /im) and inclined to the line of sight by 42.8+^5 
degree. Morever, with a very simple model, we were able to 
measure the flaring index which was found to be 1.26±%° 5 5 . 
In this model, the PAH-emitting region is only located at 
the surface of the disk, whose surface scale height H s varies 
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with v q (q is the flaring index) and whose flux intensity 
I follows r~ p . It is only recently that images of the flar- 
ing struct ure of the disk around Herbig Ae obje cts became 
available (|Perrin et al.ll2006t lLagage et al.ll2006D . It is now 
primordial to confront the spatial distribution constraints 
with existing full disk models only based on fitting the SEDs 
of the Herbig Ae stars. 

In this paper, we are interested in particular to check 
if these new observational spatial constraints can be ac- 
counted for by predictions based on existing models of pas- 
sive centrally irradiated hydrostatic disks made to fit the 
SEDs of the Herbig Ae stars. We explore whether a more so- 
phisticated disk model could explain the SED of HD 97048 
and at the same time, the spatial distribution of the cir- 
cumstellar material with little free par ameters. To the best 
of our knowledge, few tests like that (jDoucet et al.l [20061 
have been done yet for disks around pre-main sequence 
st ars of intermediate mass . To do that, we used the mode l 
of iDullemond et all (j200lD : lDullemond fe Dominikl (|2004al ) 
which could already account for the global shape of the 
SED of a quite large number of Herbig Ae stars. The dust 
model takes into account grains at thermal equilibrium and 
stochastically-heated PAH. The structure of the disk is cal- 
culated with hydrostatic equilibrium and a radiative trans- 
fer in two dimensions is used to calculate the emission of 
the different grains population. 

The paper is organized as follows. In Sect. 2, we present 
the knowledge on HD 97048 and its circumstellar material. 
In Sect. 3, we describe the observations and the data results. 
In Sect. 4, we describe the disk model and in Sect. 5, we 
compare the observed SED and spatial distribution of the 
circumstellar material to the model predictions. 

2. HD97048 

HD97048 is a nearby HAe star of spectral type Be9.5/A0 
located in t he Chamaeleon I dar k cloud , at a distance 
of 18 pc (|van den Ancker et all (|1997t ). IWhittet et all 
(|l997h ). It is surrounded by a large amount of circumstel- 
lar material left from the star formation process, which 
produces a large infrared (IR) excess over the stellar 
emission (large IR~0.35-0.40L + |Acke fc van den Anckerl 
(p00l : IVan Kerckhoven et all (j2002j)). HD 97048 has 
been classified as a HAe star o f group I with evidence 
of flared di sk dMeeus et all l200lh . s i nce it s SED is rising 
in the IR (|Acke &; van den Anckerl I2004D . Spectroscopic 
observations of the IR excess have also revealed the 
presence of strong PAH features at 3.3, 6.2, 7.7, 8.6, 11.3 
microns and nano-diamond s features in the 3.4-3.5 /xm 
regio n (|Siebenmorgen et"aTI 120001 : I Van Kerckhoven et al.l 
2002'!. No silicate emission band at 10 /im appears in the 
spectra of HD 97048. Recent mid-IR long slit spectroscopic 
observations with TIMMI2 show that the aromatic emis- 
sion features at long wavelength (i.e., 8.6 and 11.3 /im) are 
extended and come mostly from a re gion of 200-300 AU, 
likely a disk (jvan Boekel et all 120041 ) . These results have 
been confirmed r ecently by imaging d ata taken in the PAH 
band at 8.6 /mi (jLagage et al.l 1200(f ). The large extended 
pa rt of the mid- I R emi ssion seen on scales of 5 to 10 arc sec 
by iPrusti et al] (|l99l and ISiebenmorgen et all (|2000D is 
most likely due to an extended envelope of transiently 
heated very small grains and PAHs surrounding the star 
and the disk system. Using Adaptive Optics high angu- 
lar resolution (~0.1 arcsec) spectroscopic observations, 



Habart et al. (2004, 2005) were also able to spatially 
resolve the emission in the aromatic and diamond features 
around 3 /jm and found that the emission must be within 
30 AU closely r e lated to the star-disk system. Finally, 
IWeintraub et al.l (|2005f ) have recently reported near-IR 
molecular hydrogen emission for this object. 



3. VISIR observations 

3.1, Observations 

The observations were performed using the ESO mid- 
infrared instrument VISIR installed on the VLT (Paranal, 
Chili), equiped with a DRS (former Boeing) 256x256 pixels 
BIB detector. The object was observed on 2005, January 
25 th and June 17 th . It was observed in the PAH bands at 
8.6 and 11.3 /im and in the a djacent continuum. T he data 
at 8.6 /tm were already used in lLagage et al.l (2006). In this 
paper, we focus on the data obtained at 11.3 /im. A sum- 
mary of the observations is given in Table [TJ These data 
were taken with an imaging mode of VISIR which allows 
diffraction-limited image in the N band: the BURST mode. 

Under good seeing (< 0.5 arcsec in the visible), the 
images in the mid-IR are diffraction-limited even on a 8 
meter class telescope. Unfortunately, the median seeing 
experienced at Paranal is rather of 0.8 arcsec, which 
degrades the angular resolution. Indeed, for a seeing of 
0.8 arcsec in the visible, assuming that the wavelength 
dependence of the seeing follows a A -1 / 5 law, the seeing 
value at 10 /im is 0.4 arcsec, which is larger than the 
diffraction limit of 0.3 arcsec. This represents a 5 pixels 
movement on the detector with the smallest field of view 
of VISIR (0.075" /pixel). In order to reach the best spatial 
resolution with VISIR, we implemented a new imaging 
mode on bright objects, the BURST mode. The principle is 
to take short enough exposure images (< 50 ms) to freeze 
the turbulence; the coherence time of the atmosphere at 
10 /im is around 300 ms at Paranal for a good seeing. In 
order to correct for the turbulence by offline processing, 
the data are stored every 1000 elementary images in one 
nodding position for a chopping frequency of 0.25 Hz 
in the direction north/south. The nodding direction is 
perpendicular to the chopping direction with an amplitude 
of 8". After classical data reduction in mid-IR, a cube 
of 5OC0 images chopped and nodded (4 beams/image) 
is obtained. Because of the turbulence, each source on 
an image moves independently and as a result, we have 
to extract individually the 4 sources in each image (4 
quarters) of the cube and shift and add the image with the 
ones corresponding to the same quarter. Finally, we shift 
and add the four final images of the four quarters. 



3.2. Results 

In this section, we present the images of HD 97048 obtained 
in the PAH emission filter, as well as in the adjacent con- 
tinuum emission, and compare it with the standard star in 
the same filters. 

Fig. [I] shows the image of HD 97048 in the PAH fil- 
ter centered at 11.3 /tm and compared with the standard 
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Table 1. Observations of HD 97048 in the different VISIR niters in imaging BURST mode. The sensitivity were calculated 
for each night in the BURST mode. The filter are free of any strong atmospheric line contribution. 
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0.16 
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25 
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0.80 


1.74 


HD91056 


PAH2 


11.26 


0.59 


4.5 


50 


800 


25/01/2005 


0.55 


1.66 


HD85503 


Nell 


12.27 


0.18 


7.7 


16 


320 


17/06/2005 


0.80 


1.74 


HD91056 





Fig. 1. On the left, HD 97048 in the PAH2 filter (centered 
at 11.3 /im) with VISIR camera. On the right, PSF refer- 
ence (HD85503) with the same filter. The two object are 
normalized in order to see the extension. The object is more 
extended than the reference with an asymmetry east/ west 
in the emission. 



star HD85503. HD 97048 is quite extended - up to 2 arcsec 
- compared to the reference star. The full-width at half- 
maximum (FWHM) of HD 97048 is 1.2 times that of the 
standard star, close to the diffraction-limit, for almost all 
filters (Tab. [2]). In addition, comparing the emission in the 
PAH band and in the adjacent continuum (Fig. \2$, the emis- 
sion in PAH is much more extended than in the continuum. 
Fig. [3] shows the brightness spatial distribution along a cut 
in the north/south and east/ west direction in the PAH fil- 
ter and in the adjacent continuum (SIV) compared to the 
standard star in SIV filter. In the PAH filter, the disk is ex- 
tended in the direction north/south up to 320 AU (surface 
brightness of 53 mJy/" 2 ) whereas the extension in the con- 
tinuum at 10.5 /im goes to 135 AU (surface brightnesss of 
110 mJy/" 2 ). These results are in agreement with those al- 
ready found in the previous study at 8.6 fim (|Lagage et al.l 
2006). Furthermore, it is also interesting to point out that 
the continuum is extended, when comparing the object and 
the reference star in the SIV filter. Finally, by comparing 
the VISIR and ISO (Infrared Space Observatory) fluxes, 
we found that the PAH emission peak to the continuum 
is stronger in the ISO spectrum. This is certainly due to 
the surrounding nebula, also inc luded in the larger beam 
of ISO, as already suggested by Ivan Boekel et all (|2004j) . 
Based on the fluxes measured in the wavelength bands at 
8.6 and 11.3 fim, we estimate that the nebula contami- 
nation is about 40%, in agreement with that found by 
Ivan Boekel et al.l (|2004l ). 



Fig. 2. HD 97048 in the PAH2 filter (image in the middle), 
and in the adjacent continuum (SIV on the left, Nell on the 
right). The images have the same signal-to-noise so that it 
is possible to compare the extension. 

4. Disk model 

In order to anal yse the observations, we use the disk 
model described in Dullemond et al. (200l[); lDominik et al.l 
(|2003D ; lDullemond fc Dominikl (|2004aD . including the most 
up to date understanding of disk structure and physics 
around Herbig Ae/Be stars. In the following, we describe 
the disk structure and the radiation transfer used and give 
the details of the adopted dust model. 

4.1. Disk structure 

We consider disks heated by irradiation from the cen- 
tral star, in hydrostatic equilibrium in the vertical di- 
rection, with dust and gas well mixed (flared disks, 
IChiang fc Goldreichl (|1997l )). In this model, the stellar flux 
impinging with flaring angle a upon the disk is absorbed in 
the upper layers of the disk, which will reradiate half of the 
flux away from the disk and half down into its deeper lay- 
ers. The inner boundary (rim) is directly exposed to the 
stellar flux and is puffed up since it is hotter than the 
rest of the disk. The model computes how much the in- 
ner rim puffs up, and how much o f the disk behind it wi ll 
be shadowed by this puffed-up rim (|Dullemond et al.ll2001l ). 
Once the star is given, the disk structure (i.e pressure scale 
height and flaring angle) is completely defined after speci- 
fying the inner and outer radii, the surface density distri- 
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Table 2. Comparaison of the FWHM for the object and the reference star (PSF) in the different filters. We also mention 
the theorical value of the FWHM in order to show that the diffraction-limit is obtained with the BURST mode. In order 
to see the spatial extension, we note in the third and fourth columns the distance from the star at 10 a above the noise 
in the east /west direction. 
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Fig. 3. Normalized intensity profiles along a cut through 
the VISIR images of HD 97048 in the PAH band (full line) 
and in the continuum SIV (dashed line) superimposed on 
a reference star (dashed-dot line) in SIV band. Direction 
East /West is in the upper panel and above, it represents the 
direction South-North. Comparing the PAH band and the 
adjacent continuum (for a same signal-to- noise), HD 97048 
is much more extended in the PAH band. Comparing the 
PSF and the object in SIV, we see that HD 97048 is also 
extended in the continuum. 



bution (£ = £o(-R/-Ro)~ p ,with Ro a fiducial radius) and 
the dust model. The models are appropriate for disks that 
are optically thick to the stellar radiation. It is the case for 
disks around pre-main sequence stars, up to very large radii 
(e.g., 5000 AU for a disk mass M d ~ O.2M , p=l). 

4.2. Radiative transfer 

We use the 3 -dimensional Monte Carlo rad iative transfer 
code RADMC (|Dullemond fc Dominikl l2004al ). for which a 
module to treat the emission from quantum-heated PAH 
molecules has been included. This module will be de- 
scribed in detail in Dullemond et al. (in prep.), but a 
rough description has bee n given by Pontoppidan et al . 
(2006). The code RADMC (iDullemond fc Dominikl l2004at 



iPontoppidan fc Dullemondl 120051 ) solves the temperature 
structure of the disk in a Monte-Carlo way using a variant 
of the algorithm of Bjorkman & Wood (1997). This Monte- 
Carlo code also produces the source terms for scattering, in 
the isotropic-scattering approximation. 

4.3. Dust model and PAH properties 

The dust is a mixture of grains in thermal equilibrium and 
transiently heated PAHs. 

— Thermal grains: The grains are composed of graphite 
and silicate with opt ical constants from Draine (1985). 
They have a MRN (jMathis et al.lll977l ) size distribu- 
tion (n(a) oc a -3 ' 5 ) with a size between 0.01 psm and 
0.3 /im. Since no silicate emission features have been 
detected in HD97048 (see Sect. 2), we have considered 
the hypothesis of thermally decoupling the carbon and 
silicate grains (See Sect. 5.3). The optical to mid-IR 
opacity ratio of silicate is much smaller than that of car- 
bon. When the grains are thermally decoupled, the low 
temperature of the silicates produces very weak features 
and the emergent flux will be dominated by emission 
from the carbon grains. But the absence of the silicate 
feature could aslo be due to a geometrical effect, to a 
lower silicate abundance, and/or larger silicate grains in 
the inner regions. 

— Transiently heated grains (PAH): We can explain 
the observed PAH spectra, at least of the isolated HAe 
stars, with PAH abundances and qualit ative proper- 
ties si milar to those of PAHs in the ISM (|Habart et alj 
l2004al) . In the ISM, PAHs are made up of a few tens 
up to a few hundreds of carbon atoms; for reasons of 
simplicity, we take only one PAH size in our model, 
Nc — 100. The hydrogen to carbon ratio is H/C — 
Jh x (6/Nc)° 5 (case of compact symmetric PAHs, see 
lOmontj ([1986)) with fn - the hydrogenation fraction 
of the molecule. Here, we consider fn equal to 1 (i.e., 
essentially fully hydrogenated PAHs) or 0.5 (partially 
hydrogen ated PAHs). We tak e the absorption cross sec- 
tion from lLi fc Drain el (|200ll ) based on both laboratory 
data and astrophysical spectra. Those authors consider 
the bands at 3.3, 8.6, 11.3, 11.9 and 12.7 psm from vi- 
brational modes of the aromatic C-H bond; the strong 
bands at 6.2 and 7.7 fan due to vibrations of the aro- 
matic C-C bonds; and a few weak features probably 
caused by C-C bending modes at 16.4, 18.3, 21.2 and 
23.1 /im. With respect to the charge, in order to keep 



C. Doucet et al.: HD 97048 : a closer look to the disk 



5 



the model simple, we assume that all PAHs are neutral. 
Simple determination of the ratio between the photoion- 
isation rate of the grain to the electron-grain recombi- 
nation rate suggests that this is probably the case in the 
outer regions ( R > 150 AU) of a disk heated by a typical 
HAe star (see lHabart et all ((200 4b ) for more details). 
Finally, we do not take into account photo-destruction 
of PAHs in a strong FUV radiation field. 

In summary, we adopt a dust model with large ther- 
mal grains of graphite and silicate and small transiently 
heated aromatic particles. The silicate abundance in dust 
is [Si/H] = 3.10 -5 , and the total carbon abundance in dust 
is [C/H] = 2.2 10~ 4 . Of this 10% are in PAH and 90% in 
large grains. 



Hs/R 
Hp/R _ 




5. Comparison between model and observations 

In this section, we compare the model's predictions to the 
observations. In the following, we first describe the adopted 
stellar and disk parameters for HD 97048 and then discuss 
the results of the calculation and the confrontation to the 
observations. 



5.1. Stellar and disk parameters 

A Kurucz model spectrum is taken for the central star 
with T e ff — 10000 K. The stellar parameters (L* = 32L Q , 
M* = 2.5Mq) were chosen from stellar evolu tionary tracks 
by S iess et al. (2000) for an age of 3 Myr (|Lagage et al.1 
2006) . Once the stellar flux has been redenned, the resulting 
SEP is in agreement wi th the photometry extracted from 
Hillenbra nd et all (jl992t ) in the UV a nd near-IR. To c orrect 
for extinction, we used the method oflCardelli et ; al] (11 989) 
where we adop t Av=1.24 (Ivan den Ancker et~ 19981 ) and 
E(B-V)=0.36 (jDavies et al]ll991t iThe et al.lll986h . 

Concerning the disk's parameters, the surface den- 
sity is taken equal to S = So (r/Ri n )~ q with q the 
pow er law index equal to 3/2 inferred for the solar neb- 
ula (| Weidenschillingl [l977h and S = 444 q. cm~ 2 a mini- 
mum value deduced from VISIR observations (jLagage et alJ 
2006). The inner radius Ri n is at the dust (silicate) evap- 
oration radius at 0.4 AU from the central star and the 
outer radius is at 3 70 AU as deduced from observations 
(jLagage et al] |200"6[) . Finally, as the disk is vertically op- 
tically thick at the wavelength of observatio n, a minimum 
mass of 0.01 M (gaz + dust) can be derived (jLagage et al.l 
2006). We take this minimum disk mass as a first guess. 



5.2. Disk's structure 

Figure |4] shows the run with radius of the pressure (H p ) and 
photospheric (H s ) scale height. The disk is flaring with H s 
increasing with radius as R 9 / 7 (jChiang &; Goldreichl[l997l ). 
At a radius of 135 AU, H s = 51 AU. 

In a previous work, we used the VISIR PAH band im- 
age to cons train the parameter s of the disk structure around 
HD 97048 (jLagage et al.ll2006j ) . Using a very simple model, 
iLagage et al.l (|2006l ) measured a flaring index of 1.26^05 
and H S =51.3*®'7 S AU at 135 AU. Both values are very close 
to those expected from our hydrostat ic, radiative equilib- 
rium models of passive flar ed disks (jChiang fc Goldreichl 
119971: iDullemond et al.ll200lh . 



Fig. 4. Structure of the disk for the template model. The 
surface (dashed line) and pressure (full line) scale height of 
the flared disk versus the distance from the central star. 



Concerning the structure of the inner region, VISIR has 
clearly not enough spatial resolution to constrain it. The 
puffed inner rim and the shadow region lie within the cen- 
tral pixel of VISIR (1 pixel is equivalent to 13.5 AU for a 
distance of 180 pc) and the two effect^ compensate each 
other in terms of resulting mid-IR emission. Finally, we 
would like to underline that as soon as the disk is vertically 
optically thick in the mid-IR, the change of the slope of the 
surface density in the model does not have any influence on 
the structure of the disk or its mid-IR emission and cannot 
be constrained here. 



J 10 




10 



100 



1 000 



A [fir 



Fig. 5. SED calculated with the template model (dashed 
line). The disk is inclined by 43 degrees and the system 
is situated at a distance of 180 pc. A Kurucz model spec- 
trum is taken for the central star with T e ff = 10000 K. 
The luminosity is L* = 32Lq and the mass M* = 2.5M Q . 
The disk is flared with a total mass of 0.01 M Q , Ri„=0.41 
AU and R out =370 AU. Full red line shows ISO SWS 
spec trum of HD 97048. Poin ts of photometry are taken 
from Hill enbrand et al.l (119921) (open diam ond black), IRAS 
(open diamond blue). iPrusti et al] |l994) (blue crosses) and 
VISIR measurement (black triangle) 



Increasing the height of the inner rim augments the shadow 
region. 
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Fig. 6. The non PSF-convolved predicted image for our 
template model in the adjacent continuum (at 10.5 /im on 
the left) and in the PAH band at 11.3 /im (on the right). 



5.3. Spectral energy distribution 

Figure [5] shows the calculated spectrum in the 0.1-1300 /tm 
range for a star/disk system inclined by 43 degree and situ- 
ated at 180 pc. At short wavelength between 0.1 to ~l/un, 
the emission comes from the stellar photosphere. The near- 
IR emission (around 3 /im) is induced by the puffed inner 
rim. The region which emits between 5 and 8 /im corre- 
sponds to the shadow region. In this part, the disk is not in 
sight of the star itself, but it receives flux from the inner rim, 
which is sufficiently strong to keep the disk up. At larger 
radii, the flaring disk reappears from the shadow, and pro- 
duces the rise observed in the spectrum around 20-30 /im. 
Concerning the dust spectroscopic features, the PAH emis- 
sion features are clearly visible, and some of them at 3.3, 
6.2, 7.7 and 11.3 /im are very strong. Strong ratio between 
30 and 10 /im and PAH features show that the disk inter- 
cepts a large fraction of energy in the outer part and are 
indeed evidence of a flared disk geometry. For geometrically 
flat disks, PAH features are predicted to be very weak, even 
when PAH with standard properties are present on the disk 
surface (Habart et al., 2004). The infrared emission under 
the narrow PAH features is mostly due to the large thermal 
grains which are very hot in the inner regions. For the model 
in which carbon and silicate grains are thermally coupled, it 
is possible to see, for example, the strong broad feature due 
to silicate emission peaking at about 10 /im. On the other 
hand, for the model in which the grains are decoupled, one 
can see (Figf5]) that this feature almost disappears. 

In Figure we also compare the predicted SED 
to one of HD 97048 constructed with different sets of 
data. Photomet ric points from th e visib le to the mid-IR 
were t aken from iHillenbrand et al.l (|1992f ) and lPrusti et al.l 
(|1994h . To that we have added IRAS and VISIR photo- 
metric measurements. Spectroscopic observations were ob- 
tained with ISO-SWS (Short Wavelength Spectrometer) 
(|Acke &: van den~A ncker 2004T). The model reproduces cor- 
rectly the global shape of the observed SED of HD97048. 
The agreement between predicted slopes and absolute 
fluxes from the near- to mid-IR waves with observed ones 
is rather acceptable (differences <20-40%). Moreover, our 
model reproduces also well the observed intensity of the 
most commonly observed PAH features, i.e., 3.3, 6.2, 7.7 
and 11.3 /im, especially if we correct for contamination by 
the associated reflection nebula (about 40%, see Section 



3). This is in agreement with previous comparisons made 
by Habart et al. (2004) between flared disk model results 
and observed PAH emission features with ISO and ground- 
based telescopes of some thirty HAcBc stars, including 
HD 97048. 

Nevertheless, one can note that there is some mismatch 
between the predicted and observed PAH spectra; the 
strength of the 8.6 /im feature is, for example, significantly 
underestimated by the model. However, this is not surpris- 
ing considering the uncertainties on the PAH absorption 
cross section (Li & Draine 2001) and our simple hypoth- 
esis that PAH are characterized by a single charge state, 
size or hydrogenation parameter along the disk. This is un- 
likely to be the case and is briefly discussed in this section. 
Also, one can note that HD 97048 presents peculiar strong 
features that peaks at 3.43 and 3.53 /im not predicted by 
our model. Several studies have proposed attributing these 
features to sur face C-H stretching modes on nanodiamond 



parti c les (e.g., Guillois et al.l Il999t [Van Kerckhoven et all 



120021 : ISheu et all 120021 : 1 Jones et al.ll2004D . Because of the 
good match between laboratory and observed spectra, this 
identification appears convincing. 

In addition, the absence of the small grains of silicate 
features in the HD 97048 spectra could appear intriguing 
since they are the most abundant dust species in interstel- 
lar space. However, in ISO spectra, Acke & van den Ancker 
(2004) reported non detection of silicate feature for 16 ob- 
jects out of 46 showing that the absence of the silicate fea- 
ture is a common phenomenon among HAeBe stars. The 
silicate emission at 10 /im is arising from grains with a size 
of 0.1 /im thermally heated in the inner region (1 < r < 10 
AU) from direct and/or indirect irradiation by the central 
star. The absence of this emission could probably result 
from various effects concerning either the dust properties 
either the disk geometry in the inner region. Here, we have 
considered the thermal decoupling hypothesis but it could 
also be to less silicate abundance in the inner disk part. 
Decreasing globally the silicate abundance could not be a 
solution since the 20 microns feature shows the presence of 
silicates in the 50 AU < R < 100 AU parts of the disk. 
We have also tested different structures of disk "made by 
hand" and we have shown that it is possible to explain the 
abscence of the silcate emission at 10 m with a geometri- 
cal effect. But at the moment, the structure is calculated 
SELF-CONSISTENTLY and the aim of the paper is not 
to construct a structure by hand that could fit all the data 
of HD97048 but won't be physical. Furthermore, this is- 
sue could not be invistigated with the resolution of VISIR 
observations. 

Finally, concerning the longer wavelength, we find that 
the submm fl ux is too low by a factor 7 compared to ob- 
servations (|Henning &: Launhardtj [l998). This shows that 
a large reservoir of large grains (around a millimeter size) 
must exist in the outer regions of the disk. These grains 
would naturally reside close to the midplane and therefore 
do not affect the shape of the disk. They are therefore not 
within the focus of this paper, and for this reason we do 
not include them in our model. If we would have included 
them, they would only affect these long-wavelength fluxes, 
because the disk has a flaring geometry (see e.g. Dullemond 
k. Dominik 2004a, Figs. 6 and 7). 
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Fig. 7. On the left, cut in north/south and east /west (up panel) for the template model of the PAH band at 11.3 /im 
convolved with VISIR PSF (dashed line) compared to the observation (full line). On the right, the same as on the left 
with the adjacent continuum at 10.5 /im. The errors drawn are 1-a errors RMS. 



5.4. Imaging 

Figure [6] shows the modelled disk image in the PAH band 
at 11.3 /im and the adjacent continuum. These emissions 
both originate from the optically thin surface of the disk 
but the PAH emission is much more extended than the 
adjacent continuum. The continuum reaches 50 % of its 
integrated intensity at a very small radius (about 2 AU), 
while the PAH feature does so at large radii (about 100 
AU). This behaviour basically reflects the different exci- 
tation mechanisms of grains in thermal equilibrium and 
PAHs transiently heated. Only grains very near the star are 
warm enough to emit at 11 /im whereas PAHs farther away 
can be excited and emit in the 11.3 /im feature. The emis- 
sion at 11.3 /im (band + continuum), in the central part 
of the disk (< 2 AU) is dominated by the thermal emis- 
sion from silicate and carbon grains, whereas in the outer 
regions of the disk, the emission is dominated by the PAH 
feature. Morever, among the PAH features, the 11.3 /im is 
one of the most exten ded spatially. Indeed, as discussed in 
Hab art et alj (|2004bf ). the features at shorter wavelengths 
are in fact stronger in the inner part of the disk (where 
PAHs are hotter because multiphoton events occur where 
the radiation field is most intense), decreasing rapidly in the 
outer cold part. On the other hand, the features at longer 
wavelengths are more extended. 

Fig [7] shows the brightness emission profiles of the PAH 
feature and the adjacent continuum obtained by convolving 
the model with the corresponding PSF observed by VISIR. 
The model predicts a spatial distribution of the PAH emis- 
sion very similar to that observed. The predicted FWHM 
and wings extension are in fact very closed to the observed 
ones (differences <20%). Moreover, our disk model predicts, 
as in the observations, an asymetry east/west, which only 
results from the inclination of a flared disk optically thick at 
the observed wavelength. This acceptable agreement gives 
a strong support for the physics underlying in our flared 
disk model. The disk parameter that most affects the PAH 
emission is in fact the disk flaring angle, which determines 
at each radius the fraction of FUV intercepted by the disk 
surface. Lower values of the flaring could be caused by a va- 



riety of reasons , for example if the dust settle s toward the 
disk midplane (jDullemond & Dominiki l2004bh . If the disk 
is less flared, the PAH emission which directly tracks the 
illumination of the disk surface will strongly reduce in the 
outer disk region. Less flared disks will have less extended 
PAH emission features and weaker. This will be particu- 
larly true for the features at long wavelength, such as the 
11.3 /im one, which have a large contribution from the outer 
disk. In the extreme case of a fully self-shadowed disks, the 
PAH feature strenghts should decrease by orders of magni- 
tude and the spatial distribution should be similar to one 
of the adjacent continuum. 

In addition, it is remarkable to note that, as predicted 
by the model, the observed spatial extension of the 11.3 /im 
feature is much larger than that observed for PAH feature 
at t he 3.3 /im, which are extended on a scale of (several) 10 
AU ()Habart et al.ll2004al) . This has the interesting implica- 
tion that PAHs appears to be present over a large range of 
radius; in other words, PAHs can survive over a wide range 
of physical conditions. Finally, concerning the spatial ex- 
tent of the 10 /im adjacent continuum emission, the model 
predicts that it is slighlty broader than the PSF but still 
agrees with the observations. 

It must be emphasized here that there are several com- 
plications which we have neglected. The most obvious is 
that we have assumed that PAHs can be characterized by a 
single size, hydrogenation and charge state. This is unlikely 
to be the case, and one can expect variations as a function 
of radius and depth in the disk. For example, PAHs are 
likely to be more positively ionized in the inner disk regions. 
Moreover, processes such as photo-evaporation or coagula- 
tion could affect the abundance and size of PAHs. In order 
to get some insight into the specific PAH properties, one 
needs spatial information of several band strength ratios. 
We are developing this study in a forthcoming paper. 



6. Summary and conclusions 

In a former paper (jLagage et al.1 1200H ). we were able to 
constrain for the first time the flaring geometry of a disk 
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around an intermediate-mass young star HD97048. These 
results were based on a very simple model making sev- 
eral assumptions such as a surfacic PAH emission, an op- 
tically thick disk, and a power-law function for the surface 
height and the intensity. In the present paper, using a full 
radiative transfer model based on predi cted disk geome- 
tries assuming hydro static equilibrium (iDullemond et all 
l200lt iDominik et aU [20031 : iDullemond fc Dominikll2004af) . 
we could: 

— j ustify the hypothe sis made in the previous paper 
(|Lagage et all u006|) . and therefore confirm the re- 
sults, e.g., the im ages calculated with a different code 
(jPinte et al.ll2006l ). 

— show that both SED and spatial distribution of the PAH 
emission and the adjacent continuum are in very good 
agreement with the model predictions for common disk 
parameters. 

— take the general agreement between observations and 
predictions as a strong support for the physical pictures 
underlying our flared disk model. 
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